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摘 要 : 以 土壤 理化 性 质数 据 建 立 传递 函数 是 间接 获取 土壤 水 力 参 数 的 重要 手段 。 采 集 巴 拉 格 尔 河流 域 土壤 0~ 
30 cm 和 30 ~ 60 cm 土 层 102 个 样品 ,测定 粒 径 分 布 . 干 容重 ` 有机质、 饱和 导 水 率 .田间 持 水 率 等 指标 ,利用 数理 统 
计 、 空 间 搬 值 及 主 成 分 分 析 等 方法 揭示 流域 根系 层 土壤 理化 及 水 力 特 性 空间 特征 ,以 土壤 饱和 导 水 率 ( 及 ) 为 例 对 多 


种 土壤 传递 函数 模型 进行 适用 性 分 析 ,结果 表明 :(1) 


以 砂 质 壤土 为 主 , 随 根 系 层 土壤 深度 加 大 , 砂 质 壤 土 含量 减 


小 , 壤 质 砂 土 增多 ; (2) 不 同 层 饱 和 导 水 率 空间 分 布 具 有 相似 性 ,土壤 深层 大 于 表层 的 K,; (3) 土壤 根系 层 30 ~ 60 
cm 土 层 的 新 建 传递 函数 模拟 效果 较 好 ,检验 参数 、6w 和 6iwss 分 别 为 073.0 和 1.62;(4) 估计 误差 的 大 小 与 土壤 


类 型 有 关 , 新 建 PTFs 对 于 砂 质 壤土 较为 适用 。 


关键 词 : 草地 根系 层 ; 水 力 参数 ; 土壤 分 类 ; 传递 函数 ; 巴 拉 格 尔 河 


土壤 是 植被 生态 系统 的 重要 组 成 部 分 ,土壤 水 
力 特性 是 土壤 水 分 运动 溶质 运 移 、 植 被 生态 研究 
的 基础 中 。 近 30 a, 土 壤 传 递 函 数 (pedo-transfer 
functions ,简称 PTFs ) 作 为 估 测 土壤 水 力 特 性 研究 的 
重要 方法 ,受到 学 者 们 的 广泛 关注 。PTFs 是 一 
种 利用 土壤 质地 、 容 重 、 有 机 质 含量 等 土壤 理化 性 
质 来 快速 .间接 估算 土壤 水 力 参 数 的 方法 。 

内 外 学 者 通过 建立 基于 不 同 理论 和 数据 库 的 
PTFs 对 土壤 饱和 导 水 率 和 土壤 水 分 特征 曲线 进行 模 
拟 分 析 “” ,常见 模型 主要 有 Campbell ™™ Cosby ™ , 
Wosteni | Saxton" | Puckett"! , Rawls PEAU” van 
Genuchten 模型 "与 Rosetta 模 型“ ,但 这 些 PTFs 在 建 
模 地 区 以 外 区 域 不 具 普 适 性 。PTFs 构建 方法 主要 
包括 线性 回归 JERE ATE RER 
数据 处 理 分 组 分 类 方法 ”等 ,进而 适用 于 不 同 地 区 
的 土壤 参数 估算 。 学 者 们 针对 不 同 地 区 或 土壤 类 
型 ,如 黄 淮海 平原 、 沙 据 -草地 相间 地 区 外 、 沙 地 、 
WAHE” KAU ERKKA IEI EA 
当地 的 PTFs。 

HEITIR ,鉴于 鲜 有 对 于 内 蒙古 内 陆 河 草原 流 
域 根系 层 土 壤 传 递 函数 的 研究 ,本 文选 取 巴 拉 格 尔 
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河流 域 为 典型 研究 对 象 , 系统 分 析 土 壤 理 化 性 质 和 
水 力 特 征 参 数 及 空间 分 布 规律 ,利用 6 种 数理 模型 
对 土壤 饱和 导 水 率 (K,) 进 行 估算 ,建立 适用 于 巴 拉 
格 尔 河流 域 根系 层 的 土壤 饱和 导 水 率 传递 函数 , 旨 
在 为 内 蒙古 内 陆 河 草原 流域 土壤 水 力 特征 和 草地 
生态 系统 研究 提供 科学 依据 。 


1 研究 区 与 研究 方法 


1.1 研究 区 概况 

巴 拉 格 尔 河流 域 (116°21’ ~ 119°31'E,43°57’ ~ 
45°23'N) 位 于 内 蒙古 锡林郭勒 盟 西 乌 珠 穆 沁 旗 境 
内 ,流域 面积 5350 km ,海拔 高 度 885 ~ 1876 m( 图 
1) ,草地 类 型 是 以 羊 草 和 针 茅 为 主 的 典型 草原 ,为 
大 | 陆 型 半 干 旱 气 候 , 冬 季 严 寒 、 夏 季 炎 热 , 年 均 降 水 
量 为 334 mm ,蒸发 量 为 1149 mm ,平均 气温 1.2 C, 
最 高 气温 可 达 37.5 Y ,最 低 为 -38.5 ,平均 风速 为 
15 km.h ,每 年 有 28 ~ 148 d 的 风速 可 达到 125 km- 
he 
1.2 数据 来 源 与 处 理 

按照 《土壤 环境 监测 技术 规范 》(HJ/T 166- 
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1995 ) 的 布点 要 求 , 结 合流 域 的 地 形 地 貌 特征 和 土 
地 利用 /植被 覆盖 类 型 ,布设 51 个 土壤 采样 点 (图 
1) ,在 草地 主要 根系 层 深 度 0~30 cm 和 30 ~ 60 cm 
处 采集 102 点 位 共 408 个 土壤 样品 ,测定 土壤 粒 径 、 
于 容重 .有机质 .饱和 导 水 率 田间 持 水 率 给 水 度 、 
毛管 含水 率 和 化 学 成 分 等 指标 。 

土壤 采样 过 程 中 ,为 保持 土壤 样品 的 一 致 性 ， 
采样 点 先 除去 地 表 浮 土 ,再 挖 取 土 壤 剖 面 ,每 个 土 
层 中 间 位 置 平 行 采 取 环 刀 样 3 个 和 土 袋 1 个 , 环 刀 
用 于 测定 土壤 水 力 特 性 参数 , 土 袋 用 于 物理 化 学 特 
征 指标 测定 。 土 壤 饱 和 导 水 率 用 定 水 头 法 测定 ; 饱 
和 含水 率 采 用 环 刀 底部 缠 纱 布 浸 水 测定 , 阁 2 次 称 
重 变 幅 小 于 2% , 即 为 饱和 状态 , 取 均 值 作为 饱和 含 


水 率 ; 土 壤 有 机 质 通过 重 铬 酸 钾 稀 释 热 法 测定 ; 土 
款 粒 径 采 用 激光 粒 径 仪 (HELOS-OASIS 型 ) 和 土壤 
筛 (规格 为 5 mm) 测 定 。 

土壤 数据 分 析 中 ,数理 统计 运用 SPSS 22, 空 间 
克 里 金 插值 及 制图 选用 ArcGIS 10.3, 主 成 分 分 析 采 
用 DPS 数 据 处 理 软件 系统 ,检验 采用 Excel 2010. 
1.3 研究 方法 

RA ` 粉 粒 和 砂粒 含量 是 土壤 分 类 命名 和 传递 
函数 模型 建立 的 主要 输入 要 素 ,根据 国际 土壤 质地 
分 级 标准 进行 分 类 ,土壤 粒 径 分 级 标准 为 : 秋 粒 ( < 
0.002 mm) , 粉 粒 (0.002 ~ 0.05 mm) ,砂粒 (0.05~2 
mm) (K 1), 

土壤 平均 粒 径 四 的 标准 差 o, 可 作为 传递 函数 


内 蒙古 自治 区 


m, 1876 土壤 采样 层 位 示意 
i 885 
1 研究 区 地 理 位 置 及 采样 点 布设 图 
Fig. 1 Location of the study area and sampling points 
表 1 国际 制 土壤 质地 分 级 标准 
Tab. 1 International standard for soil texture classification 
质地 名 称 黏 粒 (<0.002 mm)/% 粉 砂 (0.02 ~ 0.002 mm)/% 砂粒 (2 ~ 0.02 mm)/% 
REWE 0~15 0~15 85 ~ 100 
砂 质 壤土 0~15 0~45 55 ~ 85 
壤土 0~15 30 ~ 45 40 ~ 55 
粉 砂 质 壤土 0~15 45 ~ 100 0~55 
砂 质 黏 壤土 15 ~ 25 0~30 55 ~ 85 
FRE 15~25 20~45 30~55 
HW) eR BEE. 15 ~ 25 45 ~ 85 0~ 40 
WERE 25~45 0~20 55~75 
RERE 25~45 0~45 10~55 
粉 砂 质 黏土 25 ~ 45 45 ~75 0~30 
ait 45 ~ 65 0~ 55 0~55 
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第 i 个 土壤 样本 计算 值 。 


模型 建立 的 主要 输入 要 素 ”, 以 增加 模型 的 有 效 性 : 


L 


(xa) 
d=- 100 


式 中 : d, 为 粒 径 区 间 中 间 值 ; f J1 DK E FY LE 
百 分 含量 。 

采用 Campbell 和 Cosby‘! Wosten''*"! Saxton! , 
Puckett“ 和 统计 回归 “等 土壤 传递 函数 进行 土壤 饮 
和 导 水 率 的 计算 与 适用 性 分 析 。 

PTFs 适 用 性 精度 评估 采用 平均 误差 Sw (med), 
均 方 根 误差 Sauss (md ') 和 相关 系数 r。r 越 接近 1， 
说 明 计算 值 与 实测 值 变化 趋势 越 相似 ; Sye 反映 计 
算 值 对 实测 值 总 体 偏离 程度 ,平均 误差 绝对 值 越 接 
近 0 越 好 , 取 值 为 正 说 明 高 佑 了 实测 值 ,反之 ,低估 
了 实测 值 ; rus 反映 了 计算 值 与 实际 值 的 吻合 程 
度 ,此 值 越 小 说 明 计 算 精 度 越 高 。 


EEC 
var(y)var(?) 


式 中 :NN 为 土壤 样本 个 数 ; y, 为 第 i 个 样本 值 ; y, 为 


(a) 0~30 cm 土屋 
00 


累积 百 分 含 量 /% 


(b) 30~60 cmt J 


m | 


累积 百 分 含 量 /% 
个 
© 


2 结果 与 分 析 


2.1 土壤 质地 及 空间 分 布 特征 

依据 国际 制 土壤 质地 分 级 标准 对 102 采 样 点 数 
据 进行 分 析 统 计 ,植被 根系 层 0~30 cm Fil 30 ~ 60 
cm 土壤 质地 基本 一 致 , 大 部 分 砂粒 含量 在 60% ~ 
95% , 黏 粒 含 量 低 于 5%( 图 2)。 研 究 区 土壤 在 0~ 
30 cm 以 砂 质 壤土 为 主 , 沿 河道 分 布 在 流域 中 部 , 面 
积 占 比 70% 以 上 ,其 次 为 壤 质 砂 士 ,分散 分 布 于 流 
域 边界 区 域 ,面积 占 比 达到 25% 以 上 ,不 到 2% 的 粉 
砂 质 壤土 零星 分 布 于 局 部 区 域 ;30 ~ 60 cm 土 层 以 
砂 质 壤土 为 主 , 占 50% 以 上 , 主要 分 布 于 流域 中 部 
和 南部 , 坏 质 砂 土 占 31% ,分 布 于 流域 的 东南 部 和 
西北 部 , 低 于 10% 的 少量 壤土 分 布 于 流域 的 西南 部 
( 表 2, 图 3)。 总 体 可 知 ,研究 区 土壤 类 型 以 砂 质 壤 
土 为 主 , 随 着 根系 层 土 壤 深 度 的 加 大 , 砂 质 壤土 减 
少 , 壤 质 砂 土 增多 ,土壤 类 型 增多 。 

植被 根系 层 K( 图 4) 高 值 区 主要 出 现在 流域 上 
游 及 下 游 的 柴 达 木 苏 木 南部 ,呈现 块 状 分 布 ,K 低 值 
区 主要 分 布 在 流域 边界 附近 , 即 河谷 两 侧 区 域 。 总 
体 来 看 ,0~30 cm 和 30~ 60 em 区 域 到 存在 相似 的 
空间 分 布 特征 ,土壤 深层 大 于 表层 的 瓜 。 沿 巴 拉 格 
尔 河 干支 流 河 谷地 形 低洼 位 置 出 现 玉 高 值 区 域 ,在 
流域 边界 附近 河谷 两 侧 出 现 低 值 区 域 , 上 游 的 到 值 


5 砂粒 粉 粒 Oo BR 


1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47 49 51 


采样 点 编号 


图 2 草地 根系 层 不 同 


上 壤 层 岩 性 组 成 百 分 含 量 


Fig. 2 Percentage of lithological composition of different soil layers in grass root layer 
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表 2 不 同 土壤 类 型 面积 百分比 


Tab.2 Percentage of the different soil types area /% 
土 层 /em 壤土 FEE REWE 砂 质 壤土 粉 砂 质 壤 土 粉 砂 质 黏 壤土 
0~30 0 0 25.49 72.55 0 1.96 
30 ~ 60 9.80 1.96 31.37 50.98 3.92 1.96 


(a) 0~30 cm 土 层 


~ 河流 砂粒 含量 /% 
一 - REDE mm 98.83 
砂 质 壤土 分 界 i 
CO 研究 区 界 
CO 子 流域 界 57.50 Q 20km 


图 3 根系 层 土壤 


(b) 30~60 cm 土 层 


砂粒 含量 /% 


@ 乡镇 99.75 
~ 一 河流 

— 壤土 - 砂 质 壤 土 分 界 

aes ag ORAR 48.52 
CO 研究 

C3 子 流域 界 0 ae 


SAIS [ASM 


Fig. 3 Spatial distribution of the soil types in the root layer 


(a) 0~30 cm 土 层 


人 


(b) 30~60 cm 土 层 


0.3 — 


图 4 根系 层 玉 空间 分 布 
Fig.4 Spatial distribution of K, in the root layer 


高 于 下 游 , 这 与 河水 ` 地 下 水 由 运动 的 冲 洲 积 作用 
有 关 ,与 空间 降水 分 布 不 均匀 .草地 植被 差异 .人 类 
放牧 强度 等 因素 有 关 ,这 种 分 布 特征 体现 了 流域 产 
汇流 下 垫 面 的 复杂 性 。 

0~30 cm 土 层 的 玉 与 土壤 粉 粒 含量 砂粒 含量 


和 标准 差 极 显著 相关 (a=0.01) ,与 土壤 平均 粒 径 显 
著 相关 (a=0.05) ,因此 ,初步 选 定 土壤 粉 粒 含量 、 砂 
粒 含量 .标准 差 和 平均 粒 径 作 为 0~ 30 cm 土壤 层 传 
递 函数 建立 的 输入 变量 ;30 ~ 60 cm JB AY KE 
壤 粉 粒 含量 人 砂粒 含量 .平均 粒 径 和 标准 差 极 显著 
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相关 (a=0.01) , 与 干 容 重 、 黏 粒 含 量 显著 相关 (ao= ”建立 的 自 变 量 ,最 终 确定 0~30 cm 土 层 选用 粉 粒 含 


0.05) ,因此 ,初步 选 定 土壤 粉 粒 含量 .砂粒 含量 、 黏 
粒 含 量 .平均 粒 径 a 30 ~ 60 cm 
土 层 传递 函数 建立 的 输入 变量 ( 表 3) ,进一步 结合 
主 成 分 分 析 ,选取 累计 贡 We 


量 标准 差 为 传递 函数 建立 的 输入 要 素 , 干 容重 和 
粉 粒 含量 为 30 ~ 60 em 土 层 传递 函数 建立 的 要 素 
(#4) , DI Campbell, Cosby , Wosten , Saxton , Puckett 
等 5 种 传统 土壤 传递 函数 及 新 建 统计 回归 的 PTFs 


表 3 不 与 物理 化 学 指标 统计 检验 


Tab. 3 Statistical test between K, and the physic-chemical indices 


Blom 岩 土 颗 粒 含量 /% 基础 物理 化 学 指标 统计 指标 
ROL 粉 粒 砂粒 干 容重 /(g: cm) 有 机 质 /% 平均 粒 径 (d，) 标准 差 (cy 
0~30 -0.105 -0.394" 0.378” 0.154 0.157 0.345" -0.449” 
30~60 -0.284 -0.482” 0.469” 0.311" -0.067 0.469" ~0.615" 
注 :* 和 ** 分 别 表 示 在 0.05 和 0.01 的 水 平 上 具有 显著 性 ( 双 尾 检验 )。 
表 4 KK. 的 物理 化 学 指标 主 成 分 提取 0~30 em,30 ~ 60 cm £E ,新 建 PTFs 都 较为 
Tab. 4 K, principal component extraction of 均匀 的 分 布 在 45° 线 上 及 两 侧 , 拟 合 效 果 较 好 ,其 他 
the physic-chemical indices 模型 均 呈现 不 同 程度 的 下 移 ， : 计算 值 低 于 实测 
Er 值 ,其 中 Saxton 相 比 其 他 模型 更 适用 于 本 地 区 土壤 
020 2 L12 13.00 89.06 饱和 导 水 率 的 估算 。0 ~ 30 cm 土 层 中 ,新 建 PTFs 模 
30-0 eT a 型 拟 合 效果 较 佳 ,r 为 0.5,RMSE X 1.7, ME 为 0, 相 
关 性 较 高 ,无 偏离 ,吻合 度 高 ;Saxton 次 之 ,r 为 0.44， 


进行 饱和 导 水 率 计算 ( 表 5)。 
2.3 模型 适用 性 评估 

以 平均 误差 Sue 、 均 方 根 误差 ruse 和 相关 系数 
r 作 为 检验 指标 对 6 种 土壤 传递 函数 102 个 样品 ( 表 
5) 进行 适用 性 评估 ,对 比分 析 Cosby, Saxton, 
Wosten , Puckett , Campbell 和 新 建 PTFs 的 计算 值 与 
实测 值 (图 5)。 


RMSE 为 1.62, ME 为 0.53。30~60 cm 土屋 ,新 建 
PTFs 模 型 拟 合 效果 最 佳 ,无 论 对 于 r、ME、RMSE ,新 
建 PTFs 都 明显 高 于 其 他 传统 土壤 传递 函数 ,相关 性 

高 达 0.73 ,对 于 已 有 传统 方程 均 提高 了 0.2~ 0.3,ME 
为 0, 无 偏离 ,RMSE 为 1.62 ,吻合 度 高 。 传 统 PTFs 
中 ,Saxton 和 Cosby 精度 高 于 其 他 函数 ,相关 性 分 别 
为 0.53、0.52, RMSE 为 2.1、2.2,ME 为 0.6、0.87;Cos- 


RS 土壤 饱和 导 水 率 传递 函数 表述 式 


Tab.S Transfer-functions of soil saturated hydraulic conductivity 


名 称 PTFs 数学 模型 
K, =3.456x 10° x(1.3/p)'™ x au 10 '-3.63x10 ?-0.025) 
Campbell 5 
| 2.5x102-3.63x10 7-6,88% 10 | . ie . ig 
b= e J +0.2x (OO 
Cosby K, =0.336 X o ee 
-sarera sor 1.562x ot +0.985 In C,- 1.332 这 VE 
RE: K; =1.13 x 10° xe 
Woten üg | 4C,2-0.333{C,+C,) 
Wet: K,=1.13x10°x Le ay 
Saxton K.=0.352 [F0477 + 2.8460, + (-0.022 + 0.0840, ~0.084C +0.199¢,*)/a 
dE %= 0: Xe 
Pucket 天 ,= 156.96 exp (—0.198C,) 
a 0 ~ 30 em: K,=~1.12 In C,-2.78 Ing, + 3.14 
新 建 PTFs f 


30 ~ 60 cm: K,=6.94 xp- 1.82 x In C,- 1.34 


注 : C,、C,、C WERE BPEL W E E (% ) 
水 率 , o, 为 平均 粒 径 的 标准 差 。 


, p 、C, 为 土壤 容重 (gcm ”) 和 土壤 有 机 质 含量 (gkg"),b 为 与 粒 径 有 关 的 变量 ,0. 为 土壤 饱和 含 
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(a) 0~30 cm 土 层 
10 


r=0.43 r=0.44 
ME=0.53 
RMSE=1.62 


ME=0.84 
RMSE=1.75 


天 计算 /md 


r=0.23 
ME=1.50 
RMSE=2.24 


Kgl (m í d’) 


2 
Kesew/(m-d") 
(b) 30~60 cm 土 层 


r=0.53 r=0.52 
ME=0.60 
RMSE=2.10 


ME= -0.87 
RMSE=2.20 


Kayal(m-d") 


/=0.42 r=0.40 


ME=1.66 
RMSE=2.76 


Kag (m d?) 


Kw md’) 


ME=1.46 
RMSE=2.13 


Ksex/(m-d"") 


ME=1.54 
RMSE=2.64 


Kw (md) 


ME=1.84 
RMSE=2.70 


r=0.50 
ME=0.00 
RMSE=1.70 


6 
Ksew/(mn-d") 


7=-0.33 
ME=1.65 
RMSE=3.25 


Kseq/(m-d-") 


图 5 不 同 土 层 饱和 导 水 率 计算 值 与 实测 值 对 比 


Fig. 3 Comparison between calculated and measured values of saturated hydraulic conductivity in different soil layers 


by 和 Saxton 分 布 相似 ,计算 值 均 略 低 于 实测 值 ,相对 
于 其 他 传统 PTFs 更 适用 于 饱和 导 水 率 的 计算 与 推 
求 ,但 因为 Cosby 只 考虑 了 土壤 黏 粒 和 砂粒 ,难免 会 
有 一 些 局 限 性 ,Saxton 的 精度 低 于 新 建 PTFs。 
Wosten FI Campell 总 体 低 于 45° 线 ,Puckett 偏 离 更 为 
明显 。 综 前 所 述 , 新 建 PTFs 更 适合 于 研究 区 饱和 导 
水 率 的 估算 。 

由 图 6 可 以 看 出 ,30 ~ 60 cm 土 层 的 到 的 计算 值 
与 实测 值 的 绝对 误差 小 于 0~30 cm, 佑 算 效 果 较 


好 。 两 层 土壤 误差 较 大 的 原因 不 同 ,0 ~ 30 cm 绝对 
误差 的 大 小 与 土壤 类 型 有 关系 ,新 建 PTFs 对 于 砂 质 
壤土 较为 适用 ,绝对 误差 高 值 均 出 现在 壤 质 人 砂 土 分 
布 区 域 。30 ~ 60 cm 误差 变化 原因 相对 复杂 ,对 于 
砂 质 壤 士 到 的 估算 较为 准确 ,对 于 壤 质 砂 士 .壤土 误 
差 较 大 。 局 部 地 区 误差 不 具有 规律 性 ,如 崇 达 木 苏 
木 东部 ,这 可 能 与 取样 所 处 的 地 形 、 下 垫 面条 件 有 
关 , 结 合 实地 观察 ,此 区 域 存在 局 部 砂 土 与 砂 质 壤 
EHH ,土壤 类 型 较 复 杂 ,理化 性 质 差 异 较 大 , 且 地 
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于 旱 
图 例 
(a) 0~30 cm 土 层 @ 乡镇 
N 一 “水系 绝对 误差 
人 C3 流域 边界 e <0.35 
C3 子 流域 边界 @ 0.35~0.85 
K, @ 0.85~1.55 
@ 1.55~2.50 


- 9.00 
0.44 


@ -2.50 


研 究 
(b) 30~60 cm 土 层 pa 镇 
N C 流域 边界 ”绝对 误差 
人 ~ 水 条 @ <0.35 
CS 子 流 域 边 界 @ 0.35~0.85 
K 14.28 @0.85~1.55 
= @ 1.55~2.50 
0.22 @>2.50 


图 6 不 同 土 层 饱和 导 水 率 计算 值 及 计算 误差 空间 分 布 


Fig.6 Spatial distribution of the calculation error 


形 坡 度 较 陡 ,新 建 PTFs 并 未 考虑 这 些 要 素 ,进而 导 
致 估算 误差 较 大 。 


3 ”讨论 与 结论 


3.1 讨论 

土壤 饱和 导 水 率 (K,) 虽 然 在 田间 测量 和 实验 室 
分 析 都 能 得 到 准确 数据 ,但 由 于 测定 方法 的 局 限 ， 
费时 费力 ,通过 实测 方法 获取 大 尺度 下 的 饱和 导 水 
率 数据 十 分 困难 ,因此 ,间接 计算 土壤 饱和 导 水 
率 具有 重要 意义 。 土 壤 饱 和 导 水 率 PTFs 建立 时 选 
择 干 容重 与 颗粒 作为 输入 变量 ,与 王 改 改 等 "的 
人 研究 结果 一 致 ,说 明 干 容重 与 颗粒 对 饱和 导 水 率 起 
主导 作用 。Saxton 相对 于 其 他 PTFs 模 型 更 适合 于 
半 干 旱 草原 对 于 饱和 导 水 率 的 推 求 , 与 黎明 扬 等 
的 研究 结果 相同 ,如 果 采 用 经 验 模 型 计算 时 ,可 以 
选择 Saxton 模型 。 砂 粒 含量 随 士 层次 度 的 增 大 而 减 
少 , 与 甘 林 等 ”所 得 规律 相同 。 估 算 K; 的 空间 分 布 
与 施 枫 芝 等 ”所 得 规律 一 致 ,不 同 土壤 层 饱 和 导 水 
率 分 布 具 有 相似 性 ,深层 土壤 饱和 导 水 率 大 于 淡 
层 。 对 于 K; 在 同一 个 地 区 ,最 大 值 与 最 小 值 差距 较 
大 ,这 可 能 是 部 分 采样 点 位 于 坡地 的 四 部 ,局 地 下 C 
松 物质 的 堆积 致使 表层 土壤 踊 松 ,容重 较 正常 偏 
小 ,孔隙 较 多 ,导致 饱和 导 水 率 呈 现 出 较 大 变化 。 
新 建 PTFs 与 孙 丽 等 所 得 结果 的 精度 相近 ,虽然 模 
型 可 用 ,但 精度 还 有 待 于 进一步 提高 。 

土壤 层 K; 估 算 效 果 与 取样 数据 受到 人 为 因素 
影响 ,地 形 ,坡度 .植被 覆盖 等 也 是 主要 因素 。 基 于 
Landsat 8 影像 对 巴 拉 格 尔 河流 域 进 行 NDVI 计 算 


(图 7) ,可 以 看 出 ,研究 区 植被 覆盖 度 较 高 ,敖包 却 
仍 、 浩 勒 图 高 勒 镇 及 达 青 宝 力 东部 及 南部 地 区 植被 
成 盛 , 巴 拉 嘎 尔 苏 木 、 崇 达 木 苏 木 及 浩 勒 图 高 勒 镇 
北部 植被 相对 稀少 , 巴 彦 衣 苏 木 附近 植被 覆盖 度 最 
低 ,总体 趋势 植被 茂盛 程度 由 巴 拉 格 尔 河上 游 到 下 
游 逐 渐 降 低 , 这 种 空间 分 布 与 土 层 玉 存 在 一 定 程 度 
的 相关 性 ,因此 ,可 以 考虑 基于 土壤 理化 性 质数 据 ， 
在 PTFs 中 融入 NDVI 数 据 , 进 一 步 提高 土壤 水 力 特 
性 参数 的 估算 精度 。 

廖 凯 华 等 “通过 对 点 估计 线性 回归 `. 非 线性 和 
人 工 神 经 网 络 进行 对 比 ,认为 非 线性 模拟 土壤 水 力 
特性 参数 效果 较 好 ,因此 ,将 此 方法 也 用 于 田间 持 
水 率 的 估算 ,估算 精 度 ( 表 6) 与 易 湘 生 等 ”相近 。 
田间 持 水 率 传递 函数 佑 算 效 果 优 于 饱和 导 水 率 土 
二 传递 函数 的 估算 效果 。 


EENAA 
河 


尔 
HERIK 
格 

拉 Es 
FPL AR RAR AA 


浩 勒 图 高 勒 镇 


图 7 研究 区 NDVI 空 间 分 布 
Fig.7 Spatial distribution of NDVI 
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函数 构建 与 评估 


表 6 土壤 田间 持 水 率 传递 函数 及 检验 


Tab.6 Transfer-functions of soil water-holding ratio 


检验 指标 
土 层 /cm PTFs 数学 模型 
r Owe Suse 
0~ 30 0.8 0 2.61 y=1.73C,-25.35p + 42.41 
30 ~ 60 0.72 0 4.58 y=-23.8C, -0.50,+6.7 In C, + 60.5 


3.2 结论 

(1) 巴 拉 格 尔 河 草原 流域 以 砂 质 壤 土 为 主 , 随 
着 根系 层 深 度 加 大 ,人 砂 质 壤土 含量 减 小 , 壤 质 砂 土 
增多 ,土壤 类 型 增多 ;不 同 层 饱和 导 水 率 空间 分 布 
具有 相似 性 ,土壤 深层 大 于 表层 的 K。 

(2) 植被 根系 层 深层 土壤 的 新 建 PTFs 佑 算 水 力 
特性 参数 效果 较 好 ,估计 偏差 的 大 小 与 土壤 类 型 有 
XK ,新建 PTFs 对 于 砂 质 壤土 较为 适用 ,偏差 高 值 主 
要 出 现在 壤 质 砂 土 分 布 区 域 ;水 力 特 性 参数 估算 精 
JE Saxton 和 Cosby 劣 于 新 建 PTFs, 但 优 于 其 他 模型 。 

(3) 草原 植被 NDVI 空 间 分 布 与 土壤 层 K, 存 在 
一 定 程度 的 相关 性 ,可 以 考虑 基于 土壤 理化 性 质数 
据 , 在 PTFs 传 递 函数 中 融和 信 NDVI 数 据 ,进一步 提高 
土壤 水 力 特 性 参数 的 估算 精度 。 

研究 成 果 对 于 内 蒙古 内 陆 河 草原 流域 生态 系 
统 保护 .降水 人 渗 .土壤 水 运动 和 溶质 运 移 等 方面 
研究 具有 一 定 参 考 ,可 以 为 社会 经 济 发 展 中 的 生态 

台 理 和 农 牧 业 管理 提供 科学 依据 。 


可 下 
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Construction and evaluation of pedo-transfer functions 
in the Balager River Basin 
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Abstract: It is important to understand the soil characteristics of the grass root layer for the restoration and 
protection of the grassland ecosystem. The pedo- transfer functions (PTFs) based on soil physical and chemical 
property data were the primary means to indirectly obtain soil hydraulic parameters. The Inner Mongolian Balager 
River Basin was selected as the typical study area, and 102 soil samples were collected in the depths of 0 to 60 cm 
of the grass roots. The indices of soil physical, chemical, and hydraulic characteristics were tested, including 
particle size distribution, dry density, organic matter, saturated hydraulic conductivity, and water-holding capacity. 
Kriging and principal component analysis were used for revealing the spatial pattern of the soil physical, chemical, 
and hydraulic characteristics parameters in the plateau inland grassland basin. Taking soil saturated hydraulic 
conductivity (Ks), for instance, the transfer functions were built and validated. The results showed that (1) sandy 
loam was the main soil, and sandy loam content decreased and loamy sandy increased with the soil depth increasing 
in the study area; (2) Ks spatial distribution of different soil layers was similar, and the Ks of deeper layer soil was 
bigger than that of the surface layer; (3) the simulation effect of the new transfer function was better in depths of 30- 
60 cm of soil layer, and ther, Ông, and Sgys, were 0.73, 0, and 1.62 respectively; and (4) the calculated error was 
related to the soil type, and the new PTFs of sandy loam was better than loamy sandy. 
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